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Abstract 

We study four dimensional non-projectable Horava-Lifshitz type gravity, in the 
case of an action with terms, cubic in curvature. For special choices of the free 
parameters of the model, we obtain two new analytic black hole solutions which 
exhibit the standard Schwarzschild asymptotic behavior in the large distance limit. 
The effect of cubic terms in the short range behavior of the black hole solutions is 
discussed. 
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1 Introduction 

A new higher order gravity model, which claims power-counting renormalizability and it 
serves as a UV completion of general relativity, has been formulated recently by Horava [1] . 
This scenario is based on an anisotropy between space and time coordinates, which is ex- 
pressed via the scalings t — > b z t and x — > bx, where z is a dynamical critical exponent. Note, 
that in the Horava-Lifshitz (HL) gravity the four- dimensional diffeomorphism invariance of 
general relativity is sacrificed in order to achieve power-counting renormalizability. 

It is worth noting that HL gravity has stimulated an extended research on cosmology 
and black hole solutions [2]-[25]. In addition to general relativity studies, quantum field 
theory models in flat anisotropic space-time have also been considered [26]-[38j. 

As it has already been mentioned, the HL gravity violates local Lorentz invariance 
in the UV, however it is expected that general relativity is recovered in the IR limit. 
This implies a very special renormalization group flow for the couplings of the model, 
but there is no theoretical study supporting this specific behavior. In addition, there are 
several other possible inconsistencies in HL gravity which have been discussed (see for 
example [39]- [45] and references therein). In particular, the absence of full diffeomorphism 
invariance introduces an additional scalar mode which may lead to strong coupling problems 
or instabilities, and in this way it prevents HL gravity to fully reproduce general relativity 
in the IR limit. However, these problems will not be addressed in the present paper. 

For the construction of the HL action, there has been proposed [TJ the so called " detailed 
balance principle". The main advantage of this approach is the restriction of the large 
number of arbitrary couplings that appear in the action of the model. However, a more 
general way for constructing the action would be to include all possible operators which are 
compatible with the renormalizability [2JH6]; this implies that all operators with dimension 
less or equal to six are allowed. 

Before proceeding it is important to mention that the HL gravity can be separated 
into two versions which are known as projectable and non-projectable. In the projectable 
version the lapse function N (see Eq. (12. 2p below) depends only on the time coordinate, 
while in the non-projectable version N is a function of both the space and time coordinates. 

In this paper we study spherically symmetric solutions for non-projectable Horava- 
Lifshitz gravity, in the case of an action without detailed balance, that includes all possible 
terms which are quadratic and cubic in curvature. For special choices of the free parameters 
of the model, we obtain two analytical solutions which exhibit the standard Schwarzschild 
black hole asymptotic behavior in the large distance limit. The effect of the cubic terms in 
the short range behavior of the black hole solutions is discussed. 

One should refer to some previous works towards this direction. In the case of actions 
with detailed balance, spherically symmetric solutions were found [I], but they had an 
unconventional large distance asymptotic behaviour. The correct Schwarzchild-flat asymp- 
totic behaviour can be recovered if the detailed balance action is modified in the IR by a 
term proportional to the Ricci scalar, and the cosmological constant term is set to zero 
[5]. A similar study, in the case of a non- vanishing cosmological constant, has also been 
carried out [6]. A generalization to topological black holes has been obtained in [7]. In 
addition, a systematic study of static spherically symmetric solutions of HL gravity (for 
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an action which includes only the quadratic terms in curvature) has been presented in [8] 
where the most general spherically symmetric solution for general coupling parameters was 
obtained. Finally, a general study of black hole solutions in the case of 5D Horava-Lifshitz 
Gravity, with quadratic terms in curvature, is presented in Ref. pU]. Note that in this 
work we consider for the first time the issue of spherically symmetric solutions for an action 
including terms cubic in curvature. 

The paper is organized as follows. In Sec. 2 we write down the most general action for 
HL gravity. In Sec. 3 we derive the equations of motion. In Sec. 4 we present the two static 
spherically symmetric solutions and finally Sec. 5 contains our conclusions. In Appendix [XI 
we show that the above mentioned two solutions, when the cubic terms are absent, coincide 
with solutions which are already known in the literature. Finally, in Appendix [B] we give 
a brief presentation of detail balance condition, and we discuss it in connection with the 
results of Appendix IA.1I 



2 Horava-Lifshitz Gravity 

In this section we introduce the notation for the so-called Horava gravity models in the 
case of three spatial dimensions (d = 3). These models are characterized by an anisotropy 
between space and time dimensions 

[t] = -z, [x] = -1 , (2.1) 

where z is an integer dynamical exponent. In order to derive the action of the model, it is 
useful to express the space-time metric in the Arnowitt, Deser and Misner (ADM) form 

ds 2 = -c 2 N 2 dt 2 + gij {dx l - N l dt) (dx j - N j dt) , (2.2) 

where c is the velocity of light, with dimension [c] = z — 1, and spatial components dx l /dt 
(i = 1,2,3). In addition, N and Ni are the "lapse" and "shift" functions which are used in 
general relativity in order to split the space-time dimensions, and is the spatial metric 
of signature (+,+,+). For the dimensions of "lapse" and "shift" functions we obtain 

[N] =0, [N^ = z - 1 . (2.3) 

The action of the model is constructed from a kinetic plus a potential term according to 
the equation 

S = TT^V / dtd*xy/\f\N {£ K + Cv} (2.4) 



WnGc , 

in which d (D — d + 1 = 4) is the spatial dimension and G is the Newton constant. 

In this paper the dynamical exponent z is set equal to 3. Note, that this choice for z is 
an immediate consequence of power counting renormalizibility request. In particular, the 
coupling \6ttGc in the above action has dimension [167rGc] = z — 3, hence, if z — 3 the 
HL gravity model is renormalizable, for z > 3 is superrenormalizable and for < z < 3 is 
nonrenor malizable . 

We emphasize that the main motivation for considering models of this type is the con- 
struction of a power-counting renormalizable gravity model. However, in order to achieve 
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renormalizibility, and simultaneously keep the time derivatives up to second order, we have 
to sacrifice the standard 4D diffeomorphism invariance of general relativity, which is now 
restricted to the transformation 

x* = x{x j ,t), t = t(t) . (2.5) 

The kinetic part in the above Lagrangian of Eq. (12.41) can be expressed via the extrinsic 
curvature as: 

C K = (K ij K ij -\K*), K ij = ^{-d t g ij + V i N j + V j N i }, i,j = 1,2,3 (2.6) 

which is invariant under the transformations of Eq. (12. 5p . Note, that A is a dimensionless 
running coupling which is expected to approach unity in the IR limit, so local Lorentz 
invariance is restored, as we explain later in detail. For the construction of the potential 
term we will not follow the standard detailed balance principle (for details see Sec. 7), 
but we will use the more general approach [21 06], according to which the potential term 
is constructed by including all possible renormalizable operators 0, that have dimension 
smaller than or equal to six, hence we write 

£y = £r + C R 2 + C r ;a + C AR 2 (2.7) 

where 

£r = Voa + ViaR, £/?2 = ri 2a R 2 + r]2bR %3 Rij (2.8) 

jC r s = m a R 3 + V3bRR ij Rij + ^ C R\R)R\ (2.9) 
Car* = V3dRV 2 R + ri 3e V \R jk V R jk (2.10) 

The dimensions of the various terms in the Lagrangian read 

[R] = 2, [R 2 ] = 4, [R 3 ] = [AR 2 ] = 6. (2.11) 

where the symbol A is defined as A = gijdidj (i — 1,2,3). In addition, we have used 
the notation R, Rij and RijM for the Ricci scalar, the Ricci and the Riemann tensors 
(i, j = 1, 2, 3), which correspond to the spatial 3D metric g^. Note that the term R^ kl Rijki 
does not appear in Cr2, as the Weyl tensor in three dimensions automatically vanishes. 

The first term Cr, with r]o a = 0, is necessary in order to recover general relativity, 
without a cosmological constant term[§, in the IR limit. The other terms C R 2, C R 3, C^ R 2 
include all possible quadratic and cubic terms in curvature, which become important in 
short distances when the space-time curvature becomes large. In addition, r] 2a , i]2b, V3a, V3b, 
T]3c, V3d, T]3 e are coupling constants with dimensions 

[Via] = 6, [V2a] = [V2b] = 4, [r) 3a ] = M = lV3c] = M = M = 0. (2.12) 

1 We have ignored the terms which violate parity, see also [46] . 

2 We are interested for spherically symmetric solutions which have the standard Schwarzschild asymptotic 
behavior, so we have set the cosmological constant rjQ a equal to zero 
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If this model is to make sense, it is necessary that general relativity is recovered in the 
IR limit. Although this difficult question has not been analyzed, we will assume that the 
renormalization group flow towards the IR leads the parameter A very close to the value 
one (A = 1), hence general relativity is recovered 0. Also, to obtain the Einstein-Hilbert 
action 

Seh = j dx°d 3 x^N - K 2 + R) , (2.13) 

without a cosmological constant term, we have to set rji a = c 2 , r] 0a = 0, and 

Ka = ± [~d o9ij + V, + V, (^pj } . (2.14) 

where the time-like coordinate Xq is defined ct. 



3 Equation of motion 



We are looking for spherically symmetric solution, in the case of non-projectable Horava- 
Lifshitz type gravity when the lagrangian part 

U R 2 = r] 3d RV 2 R + r ]3e V l R 3k V i W k 

vanishes, or equivalently 773^ = r] 3e = 0. As we will see, this assumption simplifies signifi- 
cantly the corresponding equations of motion. For the metric we make the ansatz: 



ds< 



-N(rydt 2 + 



2 '-'- 2 ' h ' 2 ■ rdil 2 



f(r) 



where r is the radius coordinate and 

dVl 2 = d6 2 + sin 2 



(3.1) 



(3.2) 



is the metric of a 2D maximally symmetric space,with spherical topology. In what follows 
it is convenient to perform the transformation 



/(r) = l + r 2 Z(r). 

The action of the model, after angular integration, can be put in the form 



S 



N(r),Z(r) 



dZ{r) 
dr 



+00 



dr L 



N(r),Z(r) 



dZ{r) 
dr 



(3.3) 



(3.4) 



3 As we have note in the introduction, a scalar mode remains, in the IR limit, due to the absent of local 
Lorentz invariance. This mode is responsible for instabilities and strong coupling problems. 

4 There is a more general ansatz for the metric, of the form ds 2 = —N(r) 2 dt 2 + /" 1 (r) (dr + N r (r)dt) 2 + 
r 2 dil 2 with nonzero shift N r (r) (see Ref. [H]), but we have set N r (r) = ( see Ref. [8]) to simplify the 
equations. 
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where 

dZ 



L 



dr 



N2 l R [r^j + P(Z) (rfj + iM(Z) (rf ) + 4Q(Z)) 



(3.5) 



and the coefficients i£, M(Z), P(Z) and Q(Z) are determined by the equations 

i?= -^(32^ + 12^ + 5%,) 

M(Z) = -(5407/ 3a + 180^ + 60t7 3c )Z 2 + (60r/ 2a + 20r? 26 )Z - 5 Vla , 

P{Z) = -(720 %a + 2507^ + 90 V3c )Z + (40r/ 2a + I5rj 2b ) 

Q{Z) = -(540r/ 3a + 180//36 + 60r/ 3c )Z 3 + (90r/ 2a + 30r/ 2b )Z 2 - 15^ la 

Note, that the coefficient R is a constant independent from r, while the functions M(Z), P(Z) 
and Q(Z) do not depend explicitly on r. In the case of the spherically symmetric ansatz, 
it is convenient to define five new couplings according to the relations 

p = 32r/ 3a + 127/36 + 5773c, V = 540t7 3(1 + 180i] 3b + 60r) 3c 

C = 720 V3a + 2507/36 + 90t7 3c , t = 307/ 2a + 1077 2 6, <p = 40r7 2a + 15t/ 2 6, (3.6) 



so we obtain 



M(Z) = -r/Z 2 + 2tZ - 5 
P(Z) = -(Z + V 

Q{Z) = -r]Z 3 + 3tZ 2 - 15Z (3.7) 



The corresponding equations of motion read: 

d ( dL\ dL _ d ( dL\ dL _ 

The first of the Euler-Lagrange equations yields for the function Zir) the equation 

* ( r f ) 3 + ™ ( r f ) 2 + ™w (?f ) + 4 «<^ 

If we solve the above equation we obtain a first order differential equation 



(3.9) 



dZ 

r— = H e {Z), £=1,2,3 (3.10) 

In general the algebraic equation ( 13 .9j) has three solutions Hi(Z), (I = 1,2,3), however a 
function H^{Z) is acceptable only if the following three conditions are satisfied 

• He(Z(r)) should be a real function for all the values of r in the range: < r < +oo 
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• H e (Z) ~ -3Z + 0(Z 2 ) for Z -> 0, this implies that /(r) behaves as /(r) ~ 1 - 2M 
asymptotically 

• N(r) 2 ~ /(r) for r — > +oo, if we demand the metric to recover the standard 
Schwarzschild form asymptotically 

Notice that there are solutions that violate the above three condition. For example, there 
are solutions with Anti de Sitter (AdS) or de Sitter (dS) asymptotic behavior. However 
they will not be examined in the present paper. We will now derive the equation of motion 
for the function N(r). If we set 

N{r) 2 



N(r) 



f(r) 



from the second of the Euler-Lagrange equations above we obtain 



dN{r) 
dr 



+ C{r)N{r) = 0, C(r) 



1 rf(r 3 Gi) G 2 



r 3 Gi dr 



rG 1 



where 



G ^ 3R { r %) 2 + 2P{z) i r %) +mz) ' 

G 2 = P'(Z) (r^\ + 4M'(Z) (r^\ + 4Q'(Z). 
It is convenient to change the variable in Eq. (I3.12p from r to Z, to obtain: 



dZ 



H(Z) 



3-& 



1 dGi 
~G x ~dZ 



where 



G X {Z) = 3RH(Z) 2 + 2P(Z)H(Z) + AM(Z), 
G 2 (Z) = P'{Z)H{Zf + AM'(Z)H(Z) + AQ'(Z). 



(3.11) 
(3.12) 

(3.13) 
(3.14) 

(3.15) 



(3.16) 
(3.17) 



Finally, we emphasize that the parameter A does not appear in the equations of motion, 
so the spherically symmetric solutions depend only on the parameters g, rj, r, ( and (p. 
The reason is that we are looking for static spherically symmetric solutions solutions with 
zero shift functions, hence extrinsic curvature terms do not contribute to the equations 
of motion, as they contain only time derivatives of the metric components (in particular 
Kij = for the metric of Eq. (13. ip ). In addition, the parameter A appears only in the 
extrinsic curvature part of the action. Therefore, the solutions we obtain in the following 
sections will be valid for arbitrary values of A. 



4 Analytic solutions 

Although, it is not possible to give a full analytical treatment of the most general case 
in the presence of terms both quadratic as well as cubic in curvature, we have obtained 
analytic solutions for two specific choices of the couplings of the model, which are presented 
in detail in the following subsections 14.11 and 14.21 
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4.1 First case: g = 0, £ = and = 

An analytical treatment of a third order algebraic equation (see Eq. ( 13. 91) ) is a difficult 
task. However, if we set p = £ = <p = 0, we can reduced this equation to a first order 
equation which can be handle analytically. In this case we have R = and P(Z) = 0, 
hence Eq. (13.91) can put into the form 

dZ = _Q(Z)_ 

dr M(Z) { ' ) 

where Q(Z) and M(Z) are given by the equations 

M(Z) = -r]Z 2 + 2tZ - 5 (4.2) 
Q(Z) = -7]Z 3 + 3rZ 2 - 15Z (4.3) 

Now, if we integrate Eq. ( 14.1 ft we find 

Z ( V Z 2 - 3tZ + 15) = -^f (4.4) 

where Cm is an integration constant which is related with the mass of the black hole, as 
we will see later. 

For the general case in which rj ^ 0, we have to solve the following third order algebraic 
equation 

7]Z 3 - 3rZ 2 + 15Z + -j- = (4.5) 

If we set 

5r]-r 2 2t 3 - 15r?7 C, 

rj 2 ' 2rf 2r\r 

the discriminant A 3 = q 3 + w 2 , of the cubic Eq. (14.51) . is 



^1 = ^-^^ w = "' 9 J""' -£3. ( 4 - 6 ) 



5r] — r 2 \ 3 I It 3 — Yhrr\ C M 



rj 2 J \ 2rf 2r]r 
After some algebraic computation we find 



25 (2O77 - 3r 2 ) X 2 + 2tCm(1 5t7 - 2r 2 )X + r^C 2 
l^X 2 " 

where we have set X = r 3 . The numerator of A 3 is a quadratic polynomial with discrimi- 
nant 

A 2 = -1QC 2 M (5 V -r 2 ) 3 (4.9) 
According to the behavior of A 3 we have three distinct cases: 

1. Non negative determinant A 3 in the range < r < +oo, equivalently A 3 = has no 
real root for X = r 3 or it has two negative real roots (r 2 < 5r/, or 5r] < r 2 < yr/ and 
tCm > 0). 
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2. One real positive and one real negative root (y?7 < r 2 ). 

3. Two real positive roots: (5rj < r 2 < an d tC m < 0). 

which are analyzed in the following three subsections. Note that although the construction 
of spherically symmetric solutions in these three cases appears different, the characteristic 
features of the solutions are similar, as discussed in subsection 4.1.8. 



4.1.1 Non negative discriminant A3 in the range < r < +00 



In this section we will study the case for which A 3 > in the range < r < +00. The 
discriminant A 3 (see Eq. (14. 8p ) is non- negative when A 2 < (see Eq. (14. 9p ) and the 
coefficient 20// — 3r 2 > 0, of X 2 in Eq. (14. 8p . is positive. These conditions are equivalent to 
5r] — r 2 > 0. An alternative way to achieve non- negative A 3 for X = r 3 > is the case of 
two negative roots, which corresponds to the condition 5r/ < r 2 < ^77 and r > 0. Hence, 
the solutions we present in this subsection are valid only when the parameters r\ and r obey 
the restrictions r 2 < br], or 5r] < r 2 < an d tCm > 0. 

In the case of non- negative discriminant A 3 , as we know from the theory of cubic 
algebraic equations, there is a unique real solution which can be expressed as 



Z(r) = s + t + 



(4.10) 



where 



or equivalently 



s = \ w + a/aT, t=\/w- x/Aj 



(4.11) 



2r 3 - Ibrrj C_ 



2rf 



M 



2r/r 3 ^ \ 



5i] — r 2 



2r 3 - Ibrrj C_ 



2r] 3 



M 



2i]r 3 



1/3 



(4.12) 



2r 3 - I5rri C_ 



2r] 3 



2r/r 3 \\ ^ 



while for the function f(r) we take 



5i] — r 2 



2r 3 - 15rr] C_ 



2r] 3 



M 



2r/r 3 



1/3 



(4.13) 



/( r ) = l + r 2 s + t + - 



(4.14) 



If the discriminant is positive there is one real and two complex conjugate solutions. If the discriminant 
is zero we have one triple real solution. 
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4.1.2 f V < r 2 

For y?7 < r 2 the discriminant A 3 has one positive and one negative root. Note that only 
in this case the coupling 77 can be negative. If r\ is the positive root of A 3 , we find that 
A3 < for r > ri and A3 > for < r < r\. We remind the reader that for negative A3 
(for r > ri) the cubic algebraic equation (14. 5p has three real solutions. In particular we 
have: 

(4.15) 

t)i+~ (4.16) 
77 

t)i + - (4.17) 
V 

with s and t defined through 

s = ype**, t = ^pe~^ (4.18) 
where p and 9 are given by the equations 



p= v /Z ?, = 




(4.19) 



For the definitions of q and w, see Eq. gS]) above. Note that q 2 < 0, since we examine the 
case of negative discriminant A 3 = q 3 + w 2 . 

We have checked that only the second solution Z 2 leads to the standard Schwarzschild 
asymptotic behavior, while the other two Z\ and Z 3 have an AdS or dS large r asymptotic 
behavior. So, for the construction of spherically symmetric solutions of the form f(r) = 
l+r 2 Z, we will use the function Z 2 for r > r\, and the function Z, Eq. (14. 10p of the previous 
section, for r < r\. Also, we would like to stress that the function f{r) is continuous and 
has a continuous first derivative on the matching point n, as required. On the other hand, 
if had used Z\ or Z 3 , rather than Z 2 , the function f'(r) would have been discontinuous at 

n. 

4.1.3 5r] < t 2 < f r] and rC M < 

In this case the discriminant A3 has two positive roots, which are denoted as r% and r 2 
(ri < r 2 ). If r is inside the interval (ri,r 2 ) the discriminant A3 is negative, while for r 
outside this interval the discriminant is positive. Similarly with the previous section we 
can construct spherically symmetric solutions, which are continuous with continuous first 
derivatives at the matching points r x and r 2 , by assuming that Zir) is given by Eq. ( I4.17p . 
if r\ < r < r 2 , and by Eq. f )4.10p . if < r < r 1 or r > r 2 . 



Z l = S + t+ - 
V 



-11- 



4.1.4 Determination of N(r) 2 

As discussed in the previous section, the Euler-Lagrange equations for iV yield 



dN{Z) + C(Z)N(Z) = 0, C(Z)- ' 



3 - -J- 
G 1 



+ J-^i (4.20) 



<iZ v ' v y ' v y #(Z) 
where 

H{z) ~ ~M{zy 

G 1 = AM{Z), 

Then we obtain 

~ ^(Z)-3M(Z) _ n 
- Q(Z) - °' 

where we have taken into account equations (14 .2p for M(Z) and Q(Z). Finally, with an 
appropriate choice of the integration constant, we find 



N(r) 2 = f(r) (4.21) 



4.1.5 Large distance limit 



For r — > +oo the asymptotic behavior of the function f{r) should be identified with the 
one of the Schwarzschild solution: 

r -> +oo => f(r) ~ 1 - — , (4.22) 

r 

where the Newton constant has been absorbed in the mass M (M MGn) . If we 
compare with Eq. (13. 3 j) we obtain the corresponding asymptotic behavior for the function 
Z(r) when r — > +oo 

r ->■ +oo Z(r) ~ (4.23) 

In order to find the asymptotic behavior of f(r) we can use Eq. (14.5p . For Z — >• 0, if we 
neglect the higher order terms Z 2 and Z 3 in Eq. (14. 5ft . we obtain: 

Z(r)^-^L + (\) (4.24) 



Now, if we compare with Eq. (I4.23P we find that 2M = (7^/15. 
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4.1.6 Short distance limit 



Note that for r — > the discriminant A3 becomes positive for all values of rj and r, hence 
the short distance limit of the function f(r) can be obtained if we expand Eq. ( I4.10p around 
the point r = 0. The result reads: 



/(r) ~ 1 - 



C 



M 



T 2 5f] — T 2 



Tj 



-r" + — 



+ O (r 4 ) 



(4.25) 



The above expansion implies a regular behavior for /(r) at the axes origin. It is worth 
comparing against the corresponding behavior of the solution in the case of zero cubic 
terms (77 = 0). Note that in the limit 77 — > we cannot use Eq. (14.2511 . as it blows up for 
rj = 0. However, in appendix A we have derived the expression in Eq. (1A.3I) for f(r) when 
rj = (Kehagias-Sfetsos solution), hence, if we expand Eq. (IA.3j) around the point r = 0, 
we obtain 



/(^l-J^+^ 2 + 0(r 7 / 2 ) 



(4.26) 



This behavior of f(r) has an essential difference with the one in Eq. ( 14. 2 5ft . In particular 
f'(r), in the case of Eq. ( 14.261) . blows up as l/y/r when r tends to zero, hence in the 
absence of the cubic terms the regular point at r = turns into a singularity. 



4.1.7 Zero quadratic terms: r = 

If we set t = in equations ( 14. 1011 . (14.121) . (14.131) . we obtain the simplified expression 



Z(r) 



1 1/3 



+ 



C 



M 



C 



M 



1r\r l 



2i]r 3 



+ 



C 



M 



2r/r 5 



+ 



c 



M 



2r]r 3 



1/3 



(4.27) 



which corresponds to zero quadratic terms in curvature. The above equation can be written 
in a more compact form if we set 

5 



C 



M 



30M, u 



The expression for f(r) reads: 



f(r) 



l + r 



'u 3 r 



6 + (3wM) 2 - 3uM) 



1/3 



V 



/(i^ 6 + {3uM) 2 + 3uMj 



1/3 



which has the following large distance asymptotic behavior 

f(r) 



2M 8M 3 ^ / 1 

1 + 77^7 + 



3ur 8 



11 



(4.28) 



(4.29) 



(4.30) 



It is worth comparing against the large distance asymptotic behavior of Kehagias-Sfetsos 
solution, Eq. ( ]A.7j) in appendix A , which is given by the following equation: 

/(r)Kl -^ + ^ +0 (^ (4.3!) 



Notice the difference in the next-to-leading terms of the above expansions. 



4.1.8 Numerical results for the first case 



1.0-1 
■ (\'\\ 




1 1 1 1 1 1 1 • 1 • 1 1 r 

1 2 3 4 5 6 

r 



Figure 1: f(r) versus r, for r = 1, rj = 0, 0.1, 1, 10, 50, 100 and C M = 40 (M = 4/3). 

In this section we study graphically the effects of the cubic terms, on the short distance 
behavior of the class of black hole solutions described in subsections 4.1.1 and 4.1.2. Note, 
that rj measures the effect of cubic in curvature terms, while r corresponds to the effect of 
second order in curvature terms. 

In Fig. [T]we have plotted the function f(r) for a fixed value of parameter r = 1, and 
several values of r] : r] = 0,0.1,1,10,50,100. Notice that rj = is the Kehagias-Sfetsos 
solution, i] = 0.1 corresponds to the solutions described in subsection 4.1.2, while the 
remaining values of r\ correspond to the solutions in subsection 4.1.1. We have chosen 
Cm = 40 =>- M = = 4/3, which is positive, so the solutions in subsection 4.1.3 are not 
shown in this figure. We see that for small values of rj = 0.1, 1, 10, there are deviations 
from the r] = solution only inside the horizon. However, for larger values of r), such 
as rj = 50 or rj = 100, when the effect of cubic in curvature terms becomes stronger, we 
observe important differences even outside the Schwarzschild horizon. Note that in general 
our solutions have two horizons as shown in Fig. (TJ but for the large value rj = 100 these 
two horizons disappear. In addition, as we show in the previous section, for rj ^ 0, the axes 
origin is not a naked singularity but a regular point. 



6 For r] = the solution is given by equation (|A.3|) in appendix A. 
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4.2 Second case: g = 0, r/ = and r = 

If we assume that p = rj = t = we have R = 0, hence Eq. (13. 9ft can be put into the form 



P (Z) ( r ^) 2 + 4M(Z)( r f) + 4Q(Z) = 



(4.32) 



where P(Z),Q(Z) and M(Z) are given by the equations 

P(Z) = -(Z + y?, M(Z) = -5, Q(Z) = -15Z 
If we solve ( I4.32p . we obtain 

dz ^10 ± 2^25 + IhfZ - \hC,Z 2 

dr if — £Z 

The above equation has the correct asymptotic behavior, H(Z) ~ — 3Z for Z — > 0, for 
the negative sign, so we will only analyze only this branch. If we integrate, assuming that 
( < (that is, setting ( = —\(\), we obtain: 



(4.33) 



(4.34) 



i + l<pZ + ^\t\zi + J^\c\z + 



15 <p 

iCfio 



15 j£_ 

CI 10 



1 + ^ + ^ICI^ 2 



10 



fZ - 1 



x 



xe 
If we set 



(^/l+f^z+fKl^-i) = Cm 



2 

M_ 
6 



= M 2 



6ICI 



50 



1 + 



15 <p 

iCfio 



' 15 y» 
ICI 10 



(4.35) 



(4.36) 



we recover the asymptotic behavior Z — > —2M/r 3 for r — )■ +00. For the function N(Z), 
we obtain: from Eqs. (13TT5D . (l3~T5j) and (l3~T7j) 



iV 2 (Z) 



c 



e Vi+|^+fici^-i 
(i + |^ + ||C|^) 



i + jW^ICI^ 2 + ^ICI^ + 



1 + + ^\(\z 2 - f-<p ( z + £ ) + 1 



15 y> 

la 10 



/ 16 jp 

CI 10 



10 



in which is an integration constant equal to 



C 



JV 



_3_^ 

10 ICI 



1 + 



15 f 

iCfio 



(4.37) 



' IS y 
ICI 10 



(4.38) 



if we demand iV — > 1 for Z — > 0, as it is required in order to achieve N 2 (r) ~ f(r) for large 
distances. Finally, we note that N(r) 2 blows up like 1/r 6 as r tends to zero (or Z — > —00), 
which implies a space time singularity at the axes origin. 
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4.2.1 Solution for zero quadratic terms 

If we set ip = in Eq. (14.351) we obtain 



i + -|CI^ 2 -i 



r<2 



(4.39) 



which can be solved analytically. The solution of the above equation can be expressed 
via the Lambert function Wl{x), which is defined as the real solution of the equation 

WL(x)e WL ^ = x, hence 



i + |lC|z»-i = wJ^ 



Simple calculations yield: 



Z = ± 



1/2 



wi 



and finally 



(4.40) 



(4.41) 



f(r) = 1 + r 2 Z = 1 -r 5 



1/2 



ici; \ 



^2 \ / ri2 



(4.42) 



Note that we have chosen the solution with the negative sign in order to achieve positive 
gravitational mass. For the large and short distance limits of f(r) we can use the properties: 



x -> =>- W L (x) ~ x 

x — > +oo =>- Wl(x) ~ ln(x) 



(4.43) 
(4.44) 



so we obtain for r — > +oo : 



/(r) 



10 V /2 Cm +0 t 1 



3|C| 



hence for the mass of the black hole we get 



M 



6ICI 



1/2 



c 



M 



which is assumed to be positive. Now, for r — > we find that: 



(4.45) 



(4.46) 



f(r) ~ 1 — r 2 



1/2 



ski; \ 



In 2 



/o2 



2 b , I 



(4.47) 



-16- 



The above expression yields a second derivative which blows up at the point r = 0, hence 
this point corresponds to a spacetime singularity. For the function N(r), if we set ip = 
in Eq. ( I4.37P we obtain 

"W' wtwm i 1+ v i+ 5 |c|Z2 J (4 - 48) 

Taking into account Eq. (I4.40p and the definition of Lambert function, we obtain the 
following closed form 

iV(r) 2 = f(r)N(Z(r)f = ) V ,^ NN 2 (4-49) 

2r*W L (yf-) (l + W L (%f-)) 

In the large r regime we find from the above equation that: 



N(r) 2 = Sir) (4.50) 



2 r 6 1 I r 12 

hence in the large distance limit we recover the standard asymptotic behavior N(r) 2 ~ /(?"). 
In addition, note that N(r) 2 blows up when r tends to zero [Z — > — oo) like 1/r 6 , as we see 
in Eq. f l4~49]) above. 



4.2.2 Numerical results for the second case 

In this section we have used the analytical expressions of Sec. 4.2 to illustrate graphically 
the black hole solutions. As we see in Fig. |2]the black hole solution has the familiar form 
with two horizons, but in contrast with the solution presented in Sec. 4.1, the function 
N(r) 2 differs significantly from f(r) inside the outer horizon and blows up for r = 0, which 
implies a spacetime singularity. Note that the solution for the first case, presented in Sec. 
4.1, are regular in the axes origin. 

In Fig. [3]we have plotted the function f(r) for a fixed value of parameter <fi = — 1, and 
several values of \(\ = 0,0.1, 1,5, 10, 100, for M — 1. We observe that when the effect of 
cubic terms becomes stronger (the parameter | £ | becomes larger) , the two horizons coalesce 
into one and finally disappear, hence for large \(\ our solution represents a naked singularity. 



5 Conclusions 



We considered four-dimensional non-projectable Horava-Lifshitz type gravity, without 
detailed balance, for an action with all possible terms quadratic and cubic in the curvature. 
Black hole solutions in the case of quadratic terms have been studied previously in the 
literature, however in this work we examine for the first time the cubic terms. Although 
in the general case it is not possible to find analytic solutions, for special choices of the 
free parameters of the model, we obtained two new black hole solutions which exhibit the 
standard Schwarzschild asymptotic behavior in the large distance limit. 
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Figure 2: f(r) and N(r) 2 versus r, for = —1, \(\ = 0.5 and M = 1. 



It is worth mentioning that the two solutions we have obtained, reduce to known so- 
lutions of Horava-Lifshitz gravity when the cubic terms are turned off. In particular, the 
solution of Sec. 14.11 corresponds to Kehagias-Sfetsos (KS) solution (see Appendix A) and 
the solution in Sec. 14.21 corresponds to Kiritsis-Kofmas (KK) solution (see Appendix B). 
In Figs. [1] and |3] the effects of cubic terms are illustrated, so the reader can compare with 
KS and KK black hole solutions in a straightforward fashion. 

In addition, we stress that the cubic black hole solutions we obtained in Sec. 14. II have 
a regular center, in contrast with the KK solution, in which the center is singular. On the 
other hand, in the case of the second solution described in Sec J4.2[ the cubic terms are not 
enough to cure the singularity (see Fig. [2]). Note, that a similar singularity is present in 
the KK black hole solution in which only quadratic terms in curvature present. 

Finally, in Sec. I4.1.7l and l4.2.1l we present the solutions when quadratic terms are turned 
off. In this case the analytical formulae for the solutions are considerably simplified. In 
addition note, that the short and large distance behavior of these purely cubic solutions 
differs from the corresponding behavior of KS and KK black holes. Further investigation 
on these black hole solutions may be a topic in a future publication. 
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Figure 3: f(r) versus r, for <p 



1, |C| = 0, 0.1, 1, 5, 10, 100 and M = 1. 



A Appendix: Zero cubic terms (g = 0, rj = and ( = 0) 

Of special interest is the case where one has ho cubic terms at all. These solutions has 
already been studied by A.Kehagias and K.Sfetsos. in Ref. [5] and by E. Kiritsis and G. 
Kofinas in [8]. 

A.l Kehagias-Sfetsos solution: (p = and r ^ 

If we set ip = in (14.41) we obtain: 



3tZ 2 - 15Z 



.3 







(A.l) 



with the two solutions 




(A.2) 



For the function / = 1 + r 



2 Z we get: 




(A.3) 
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where we have kept the negative sign to have the standard 1 — 2M/r behavior asymptoti- 
cally. Also, note that tCm should be positive, otherwise the solution has no short- distance 
limit. 

If r — > +00 we find the following asymptotic formula 

hence for the mass of the black hole we obtain 

M = %. (A.5) 



30 

This relation, along with the notation 



2t 

transforms ( 1A.3I) into: 



00 = ^-, (A.6) 



f{r) = l + ur 2 -^/r {u 2 r 3 + 4uM), (A.7) 
which is the standard form as presented in [5]. For the function N(r) we know that 

N(r) 2 = f(r) (A.8) 

Note that the recovering of the SK solution for (p = is expected, as the action we consider 
in Sec. 2 when the cubic terms are absent, is identical with the action for the detailed 
balance condition for A = 1. In the following Appendix |B] we give a detail analysis on this 
topic. 

A. 2 Kiritsis-Kofinas solution: r = and cp ^ 

If we set C = in Eq. (14. 34ft we find 

dZ 1 

■r 



,dZ = l / 1Q ± 2 / 25 - 15 z \ , A g s 

dr lp \ / 



If we integrate the above equation, in the case of the negative sign, which corresponds to 
asymptotically flat space-time, we obtain: 

yJl + l<pZ-lj ev^-^ (A.10) 

where Cm is a positive constant of integration, related to the mass M of the black hole, as 
we will see. From the definition of Lambert functions, as the real solution of the equation 
W L {x)e w ^ = x, we find 



l + l V Z-l = W L (^\, (A.ll) 



5 \ r 
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which yields: 

Z=-[\Yi[^\ +2\Y L [ -M ). (A.I2) 



while for f(r) we get 



5r 2 /„,o / Cm\ / Ca 



'< r > = 1 + ^ ^2 (^J + ^ (7? ) I (A. 13) 

In order to obtain the large and short distance limits of f(r) we can use that: 

x ->■ Wz(a?) ~ a; (A. 14) 

x -»■ +00 =>- W L (x) ~ ln(x). (A. 15) 

For r — >■ +00 we obtain 

/(r)=!l+ »^ + (J_) (A,6) 
and the mass of the black hole reads: 

M = (A.17) 
which is assumed to be positive, or ip < 0. For r — > we find that: 

^f ln2 f%y +21n (£« 

Sip \ \ r A J \ r 



/M* l + £r I In' + 21n ^) I (A.18) 



The above expression yields a second derivative which blows up at the point r = 0, hence 
this point corresponds to a spacetime singularity. For the function N(r), if we use Eqs. 
( gig) , Q33BD and f lXTTj) we obtain 



ev/i+M-i 

AT(Z) = ; (A.19) 

If we take into account Eq. (lA.lOj) and the definition of the Lambert function, we obtain 
the closed form 

iV(r) 2 = /(r)iV(Z(r)) 2 = — °\ f{r) (A.20) 

In the large r regime we find, from the above equation, that: 

N(rf =/(r) (l-2^ + O (^) ) (A.21) 

hence in the large distance limit we recover the standard asymptotic behavior N(r) 2 ~ /(r). 
In addition, note that N(r) 2 blows up when r tends to zero like 1/r 6 , as we see in Eq. ( 1A.20j) 
above. 
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B Appendix B: Comparing with the detailed balance action 

In this appendix we present briefly the "detailed balance condition", which a method 
for the construction of the potential term in the Lagrangian of HL gravity, see Sec. 2. It is 
inspired by nonequilibrium statistical physics, and it was proposed in Ref. pQ. The main 
advantage of this approach is the restriction of the large number of arbitrary couplings that 
appear in the bare action of the model. In addition, if the action of the model has been 
constructed by using the detailed balance, from an original d dimensional action, then the 
corresponding d + 1 field theory model has improved renormalization properties. These 
properties have been inherited from the original d dimensional action, see for example Ref. 
[UJ and references there in. 

In the case of the detailed balance condition the potential term, of the d+1 dimensional 
action, is derived from a variation principle, or equivalently by taking the square of the 
equation of motion of an original Euclidian d dimensional action. Particularly, in the case 
HL gravity model the action reads: 

K 2 



S v = — J dtd d x^\g\NC v , d = 3 (B.l) 

in which k 2 is a coupling constant, where 

C v = E^g irM E kl , E*= WV 9kl] ( R 2) 

and Qij-ki is the inverse generalized De Witt metric which is defined as 

giW = 1 ( g ik gj l + g il gj lj _ Xg i jg kl (B 3) 

where A is a new dynamical coupling, underlie to quantum corrections. The standard choice 
for the three dimensional (d=3) action W^/] is 

W[g M ] = W c [g kl ] + W EH [g kl ], W c = -L j d 3 x u 3 (T), W EH = fi J d 3 x^/\g~\(R - 2A W ) 

(B.4) 

in which 

and w 2 is a dimensionless coupling, while \x and Ky/ are couplings with dimensions \p] = 1 
and [A w ] = 2. Note, that the Cotton tensor, which is defined as 

Ga = e lkl V k (r{ - ]rs{) , (B.6) 



4 

can be obtain from the action Wc more specifically 

1 SW c [g kl ] 1 



Sgu w 



2 C a ( B - 7 ) 
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For the tensor in Eq. (1B.2j) . we find 

Ea = - | (& - + Aw9 ij ) (B.8) 

Now, for the potential part of the Lagrangian (see Eqs. (IB.ip . HB.2[) and (1B.3[) ) we obtain 



C v = -—djC* 3 + ^e ljk RiiVjR l k - ^-RijR lJ + 4 ^_ 3 ^ ( — ±— r2 + AwR ~ 3A w 

(B.9) 

As we are looking for spherically symmetric solutions without a cosmological constant (the 
spacetime should be Minkowsky asymptotically) it is reasonable to set A w = in the above 
action. Moreover in order to recover general relativity in the IR limit we have to add in 
the action a relevant term proportional to fi 4 R. This term violates the detail balanced in 
general, but this term becomes unimportant in the UV where the theory still satisfies the 
detailed balance condition. Finally, we take 

A< = —c^' + - ^R-' + ! §^§ j R 2 + ^ (aw) 

The corresponding full action with the kinetic term can be constructed as it is presented 
in Sec. 2. Comparing with Eq. ( 12. 131) in Sec. 2, which is valid in the IR limit, we take: 



A = l, c 2 = ^, G=£- (B.li; 

6Z71C 



We stress that in the case of spherically symmetric solutions, more specifically for a metric 
of the form of Eq. (13. II) . or 

ds i = - N {rfdt 2 + ^-+ rW (B.12) 
f( r ) 

the first two terms in Eq. ( IB. 101) . which determine the UV behavior of the model, vanish 
identically, or we have 

- — A C l3 C^ = 0, JL £ iok RN ^ = o (B.13) 

Now, we can compare straightforwardly, the Lagrangian of Eq. (IB. 101) without the terms 
which are proportional to C^C 1 ^ and e^ k RuVjR l k , with the Lagrangian of Sec. 2 that reads: 

C v = Vla R + V2a R 2 + V2b R ij R i3 (B. 14) 

when the cubic and derivative terms are absent (£r3 = and Car 2 = 0, see Sec. 2), then 
we find 

^ 2 /A1-4A) 
v* = - T , m a = 16(1 _ 3A) (B.15) 

For A = 1, we obtain that 

r}2a _8 
r}2b 3' 
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or equivalently the coupling ip = 40r] 2a + 15r] 2 b is equal to zero. We conclude that the 
condition ip = corresponds [j to the detail balance condition with A = 1, therefore the 
solution we found in Appendix IA.lt i n the case of zero cubic terms, is expected to be 
identical with SK solution of Ref. [5]. 
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